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The coherent propagation of elastic waves in a solid filled with a random distribution of pinned 
dislocation segments is studied to all orders in perturbation theory. It is shown that, within the 
independent scattering approximation, the perturbation series that generates the mass operator 
is a geometric series that can thus be formally summed. A divergent quantity is shown to be 
renormalizable to zero at low frequencies. At higher frequencies said quantity can be expressed in 
terms of a cut-off with dimensions of length, related to the dislocation length, and physical quantities 
can be computed in terms of two parameters, to be determined by experiment. The approach used 
in this problem is compared and contrasted with the scattering of de Broglie waves by delta-function 
potentials as described by the Schrodinger equation. 


I. INTRODUCTION 


In recent years, the interaction of an elastic wave with 
a dislocation has been studied within the context of con¬ 
tinuum elasticitjli^ (hereafter I). These results have been 
further used to obtain results, using multiple scattering 
theory, for the coherent propagation of elastic waves in 
the presence of many, randomly distributed, dislocation^ 
(hereafter II). Leading order perturbation theory was 
used to obtain formulae relating change in the speed 
of wave propagation to dislocation density. These for¬ 
mulae have been used to show that ultrasound—more 
specifically resonant ultrasound spectroscopy (RUS)— 
can be used as a quantitative probe of dislocation density 
in aluminurrP, with comparative advantages over X-ray 
diffraction (XRD) and transmission electron microscopy 
(TEM)P. 

The advent of a new, quantitative, nonintrusive, tool to 
probe dislocation density raises hope for progress in solv¬ 
ing long-standing challenges involving the plastic behav¬ 
ior of materials, such as fatigu^, or the brittle-to-ductil^ 
transition. More specifically, there are a number of ap¬ 
plications that stand to benefit from a nonintrusive mea¬ 
surement of dislocation density. For example the plastic 
deformation in torsion of ice single crystals has been stud¬ 
ied using hard X-ray diffractiorP in order to ascertain the 
role of geometrically necesary dislocations, and size ef¬ 
fects have been unraveled through creep measurement^. 
Also, high resolution extensometry experiments carried 
out on copper single crystals in tension have uncovered a 
rich spatiotemporal structur^. Distinct scales of plastic 
processes have been found in austenitic FeMnC steel, also 
using extensometrjffSl Finally, the Portevin-Le Chate- 
lier effecfpL -in which plastic instabilities generate local¬ 
ized bands—seriously hampers the use of some alloyJi^ 
and remains, to a large extent, a mystery. Recent at¬ 
tempts at developing a conceptual framewok to under¬ 
stand the formation of patterns by dislocations include 
those of Sethna and collaboratorsiSl, Limkumnerd and 
van der GiesserP^, and Rickman, Haataja and Le SaP^. 


Ultrasonic waves penetrate deep into a material and 
appear to be the ideal tool to probe the effects described 
in the previous paragraph. However, in order to probe 
structure, the probing wavelength must be comparable to 
the structure length scale. The results reported by Mu- 
jica et al.l^ are an average over a whole sample, and rely 
on leading order res ults in a long wavelength approxi¬ 
mation of the theorjI^H. Although these leading order 
results, obtained using perturbation theory in a multi¬ 
ple scattering framework, do provide precise and useful 
formulas, higher order approximations will be needed to 
probe shorter length scales. But, as discussed below, 
higher order results can diverge at high frequency be¬ 
cause of the zero thickness of the strings used to model 
dislocations. This paper is devoted to the elucidation of 
this state of affairs. 


II. INTERACTION OF ELASTIC WAVES WITH 
DISLOCATION SEGMENTS. 


We shall use the notation of I and II: An homogeneous 
and isotropic, linearly elastic, medium of mass density p 
is described by displacements u{x, t) away from an equi¬ 
librium position X at time t. In the absence of disloca¬ 
tions the dynamics is governed by the wave equation 


d'^Ui d'^Uk _ 


( 1 ) 


where the elastic constants tensor is given by Cijki = 
^SijSki + pt{SikSji + SiiSjk) with (A,/J,) the Lame con¬ 
stants. Longitudinal (L) and transverse (T) waves prop¬ 
agate with speeds Cl = a/(A -I- 2p)/p and Ct = a/ p/p, 
respectively. We shall call their ratio ^ = cl/ct > I- 
A pinned dislocation segment is described by its po¬ 
sition X(s,t) as a function of a Lagrangian parameter s 
and time, with orientation provided by the unit tangent 
f = X'/\X'\, where a prime denotes derivation with re¬ 
spect to s. We shall consider unbiased edge dislocations 
of length L and Burgers vector 6, that is, their position 
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at equilibrium are straight lines. In the absence of exter¬ 
nal loading the dislocation dynamics is given by a linear 
string model 

mXkis,t) + BXkis,t) -TX”is,t) =0 ( 2 ) 


where overdots denote derivation with respect to time, 
and the associated boundary conditions of pinned ends 
are Xfc(±L/2,t) = 0. In Eqn. ([^ the coefficient 

TO = ^(1-H7"'^)ln((5/5o), (3) 

defines a mass per unit length (with S and Sg the long- 
and short-distance cut-off lengths, respectively), 

r^^(l-7-2)ln(J/5o) (4) 

is a line tension, and i? is a phenomenological viscous 
drag coefficient. 

When a wave propagates in the presence of disloca¬ 
tions, there is an interaction. The behavior of the wave 
in this case is given, in the frequency domain, by the 
following equatiorPfor the velocity Vi{x,uj): 


- puj'^Vi{x,uj) - Cijkl 




dxjdxi 


Vkix,u>) = VikVk{x,uj) (5) 


where 


V^k = A M,j—6{x-X^) Mzfc — 


‘ dxi 


( 6 ) 


x=Xo 


A = 

VAij = 


8 {phf S{uj) ^ 

2 2 ^ 
TT^ m 

"tiTlj 'tjTl'i ^ 


(7) 

( 8 ) 


with n = f X t and i = b/\b\ is the unit Burgers vector 
that indicates the direction of glide. In addition, 


Sjoj) ^ y. _1_ 

iV odd ^ 


(9) 


where 


UN = 


Ntt 

~L~ 



( 10 ) 


are the eigenfrequencies of a vibrating string of length L, 
mass density m and line tension T. Only glide motion, 
that is, along £, is allowed. The potential (|^ is a sim¬ 
plified expression that captures the dislocation dynamics 
for frequencies smaller than the first resonance frequency 
wi, thus accordingly we shall use the approximation 


S{uj) _ 1 

Lu^ — uf + ioj B / m) 


( 11 ) 


When a single dislocation is present (see Fig. [^, stan¬ 
dard scattering theory can be used to turn Eqn. ^ 


into an integral equation that involves the Green func¬ 
tion (impulse response function) for the medium under 
consideration and it can be solved, for example using a 
Born approximation scheme. This has been done previ¬ 
ously for a whole spac^ and for a half space with a free 
surfac^i^. 




tr- 



FIG. 1: The basic scattering mechanism of an elastic wave 
by one oscillating line dislocation: an elastic wave of wave 
vector ki is incident upon a dislocation segment of length L 
that oscillates in response. As it does, it re-emits waves with 
scattered wave vector ks. 



A. Perturbation approach to multiple scattering 



FIG. 2: Configuration for the study of multiple scattering: the 
incident wave has transverse (T) and longitudinal (L) polar¬ 
izations and propagates through an ensemble of randomly dis¬ 
tributed and oriented dislocation segments with pinned ends. 
Each segment is of variable length L. Reproduced from Re# 

When many dislocations randomly distributed are 
present, as in Figure one studies the average (over all 
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possible configurations with prescribed probability distri¬ 
butions) displacement of the elastic medium rather than 
the displacement itself. The problem to be solved is to 
find effective plane wave solutions, with a dispersion rela¬ 
tion between frequency and effective wave number. Said 
dispersion relation is obtained from the poles of (G) (the 
brackets denote ensemble averaging) in Fourier space, 
where G is the impulse response (Green’s) function of 
the dislocation-filled medium, the solution of 


III. SCATTERING BY A SINGLE 
DISLOCATION SEGMENT: BORN SERIES TO 
ALL ORDERS 

As usual in scattering theory, we consider Eq.([^ and 
Eq.§ as an integral equation 

v^{x) = v°{x) + J dx'G^jix - x')[VjiVi]{x') (17) 


puj^G imiX-i^') Cijkl 




dxjdxi 


Gkmix^ ca) — 


^ ^ ^ikG kmi^X^ Uj') Simd(x') . (f^) 


disloc. lines 


To this end the average Green’s function is written as 
(“Dyson equation”) 

(G)-i = (G°)-i - E (13) 

where G° is the Green’s function of the medium without 
dislocations (i.e., “free”). For an infinite, homogeneous 
and isotropic medium it is, in Fourier space. 


with v^/.(x), for instance, an incident plane wave and G^ 
the free space Green function whose Fourier transform is 
given in Eq.( 16) and Eq.( 15). The T—matrix determines 


the full scattering problem 


because it satishes 


Ui(x) = v°(x) + J dx'G^j (x - ^)[Tjiv'^]{x'). (18) 


We introduce the dehnition of the T-matrix in momen¬ 
tum space through 


T,fc(fc,fc')= j dxdx'e-^^'^Tiu{x,x')e^^''^' (19) 

to obtain 

T,kik,k') = Tg\k,k')+T^l\k,k')+T^l\k,k')... (20) 


[(G°) + pk^)5ik ++p)kikk (14) 


or 


Gl{kM = 


pc^ik"^ — k^) 


{5ik—kikk)-\- 


pclik^ - kl) 


ki kk . 


, (15) 

with an implied small imaginary part in the denominator 
to insure causality. It will be explicitely considered later 
on. The corresponding expression in coordinate space 
is defined by 


G°t.(f - f') = 


(27r)3 


dfce-*'"’"G)'fe(fc)e 


ik-x 


(16) 


and we shall omit the frequency argument w in what 
follows. In the independent and weak scattering 
approximatiorfi^ E = n{T), with T the scattering T- 
matrix and (•) an average over the internal variables of 
the dislocation segment: Length, orientation, and Burg¬ 
ers vector. Their position has been assumed to be uni¬ 
formly distributed with density (number per unit vol¬ 
ume) n. IiP we investigated the properties of elastic 
waves propagating in this disordered medium and ob¬ 
tained expressions for the effective velocities and atten¬ 
uation coefficients. These quantities were obtained from 
E = n{T) at a leading order Born approximation for 
T. We now revisit this series expansion, locate divergent 
terms, and show that the full Born series for T can be 
summed. 


with the first Born approximation 

= -AM,jkjk'iMik. ( 21 ) 

The second Born approximation is computed in Ap¬ 
pendix The result is 

Tl^\k, k') = A^M,jk,Ik'iMik. (22) 

where 

I = ^ I dqq^M ■ G°(q) • M)g (23) 

Stt 1 IbTT 1 

5 pGf (g2 - 15 pGliq'^ - kl) 

with g ‘ the transpose of q. The second line is obtained 
using the properties of M: 

Mmm = 0, MifcMife = MfeiMifc = 2. (24) 

Note that / diverges because it is the integral of a 
quantity that grows ~ q^ at high wavenumbers. This 
is a consequence of the fact that we use a continuum 
mechanics approach, in which there is no intrinsic cut-off 
length. 

The third order Born approximation to the T matrix 
is worked out in Appendix]^ The result is: 

Tlf (fc, k') = -A^M,,kjI^k'iMik (25) 
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It is apparent from (251 that a power series structure isotropic on average we hav^ 


suggests itself. To pursue this possibility we consider 
the p-th term in the perturbation series and write it, 
symbolically, as 


= {VG^VG^ ■■■¥) 


ik 


(26) 


where the potential V appears p times. Replacing G° by 
its Fourier transform and repeatedly using (24) as well 


as the properties of the delta function, leads to 
T^P^k, k') = -AM,jk, (-AI)P-^ k'lMik , 
a power series, as was hoped, and formally we have 


— ^ kpl^pk 

(28) 

1+AI’ 

(29) 


With given by (21). From the T matrix it is possible 
to obtain, following the standard definitions, scattering 
amplitudes and scattering cross sections. For instance, 
the longitudinal to longitudinal scattering amplitude is 


IlUx) = (a 




-A 
1 + AI 


{ko ■ M • ko) (30) 


where x is the scattering direction and fcg the incident 
direction. Eqn. (30) is a simple generalization, to all 
orders, of a similar expression, obtained in I to first or¬ 
der. The only difference is the new denominator. The 
scattering amplitudes for mode conversion, as well as for 
transverse to transverse scattering, follow similarly, as 
does the total scattering cross section using the optical 
theorem. 


IV. MASS OPERATOR 

We return now to the evaluation of the mass opera¬ 
tor E = xi{T), where the average is over the internal 
variables of the dislocation segment: Length, orientation, 
and Burgers vector. Following II, we consider edge dislo¬ 
cations, all with the same length L, with a Burgers vector 
of hxed magnitude but randomly oriented. Dislocation 
position has been assumed to be uniformly distributed 
with density (number per unit volume) n. Thus 


^ _ -nAMij 
“ (1 + AI) 


(31) 


where 


and the integration over C is the average over internal 
variables. If we take a distribution of dislocations that is 


{Mij) = — (fcifcj + 3(5y)fc^ 

= - kkj)k'^ + ^kikj 


(33) 


where the brackets denote angular average. There are 
two objects of interest: One is the mass operator (31) 


itself. Without the density factor n and without the av¬ 


erage indicated in (32) it corresponds to the scattering 
amplitude by a single dislocation segment, the T-matrix. 
Another is the Dyson equation (13), where the mass oper¬ 
ator appears, and which provides, through the vanishing 
of a determinant, an implicit relation between frequency 
and wavenumber, from which (dispersive) effective phase 
velocities and attenuations can be read. 

Consider then the quantity 




A 


l + AI 


(34) 


that appears in the mass operator (31). The object of 
interest is the quantity / defined by (23). Introducing 
the change of variables 


Dj' = ctP 
= CLP, 


and using 


lim- = V —h iirS^x) 

X + ip X 

where V means Cauchy principal value, one has 
1 247r f 


(35) 


(36) 


I = 


(27r)3pc^ 15 

1 167r 

”*"(27r)3pc| 15 


Dp — Lo'^ — ip 


driL^i 


Dp — ui^ — ip 


1 




with 


IStt^ 

= In + 'ill 


= lo^ 
h = /oS 




1 

pcr 


dnn* 


D^ — — ip 


(37) 


dnn^ 


^n 


= hV 


D^ — — ip 

(iDD"^ uAi 

— u)'^ — ip ° 2 


dDD^ 


and 


(38) 


(39) 


(40) 


We see that I-ji, the real part of /, diverges since its 
integrand grows as D^ at large D. The imaginary part 


r 


. 1 

[37^ + 21 

/ dCMipMqjkpkq 

(32) 

° “ 157r2 

ry5 
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Ix is finite. Consider the divergent quantity I-jz, given by 
((Ml). Let 


A. Introduction of a high frequency cut-off 


In = lim /7 ?,(w,A) 

A—>-oo 


with 


IniLo,A) = Ior [ 
Jo 




iQ — u!^ 

This integral can be calculated; it is 

/ n , w 

Ini^, A) = Jo I Aw + —— w Arctanh— 

y O I\. 


for A > w. 

Calling In{<^i) = Ini this leads , in the limit A 
to 


(41) 


(42) 


(43) 


Ini — In 2 In 2 — Ino 


= 0 . 


UJT — Wn 


w; — Wo 


(44) 


So the difference between two weighted differences of I-ji, 
as specified by (44), is finite, indeed it is zero, although 


Iti diverges. It is not unusual in physics to have a finite 
effect appear as the difference of two infinite quantities; 
the Casimir effect?^ is a notorious example. So, it is 
natural to interpret ( [44| as saying that it provides the 
value of Ixi at a frequency W 3 in terms of the value of 
Iti at two other frequencies wi and W 2 . These two values 
must be independently established. 

The form of the interaction in Eqns. (5][6) suggests 


that the Born series is a long wavelengths, or low fre¬ 
quency, approximation. Indeed, since the interaction po¬ 
tential as two gradients, high wave numbers will be more 
strongly affected than low wave numbers. We should 
then have 


lim /(w) = 0 


(45) 


a condition that is already satisfied by the imaginary part 
Ix- So, suppose we wish to have In{<-^) at a finite fre¬ 
quency; consider wi,W 2 small enough that I vanishes. 
Eqn. (44) provides 


In{^) — 0 , 


in the sense that In{^) is negligible with respect to /i(w) 
at low frequencies. This last point will be elucidated in 
the next section, by a proper scaling through the intro¬ 
duction of cut-off regularization (Eq. 48). Then, from 


Eqs.(34|4'6) we have, at low frequencies, 


0 _ 

^low UJ 


A 


I -f iAIx 


(47) 


This is an expression that provides a (finite) correction, 
valid at low frequencies, to the first order Born approx¬ 
imation UBorn ~ la order to be more precise about 
what we mean by “low frequencies” we now introduce a 
cut-off dependent regularization. 


The reason I-ji , given by (38) diverges is that it consid¬ 


ers all frequencies, up to infinitely high frequencies, which 
means all wavelengths, down to infinitely short wave¬ 
lengths. This is an artifact of the perturbation theory 
approach, in which each term in the perturbation series 
introduces an integration over an infinite frequency inter¬ 
val. More precisely, the order n term has n — 1 integrals 
that extend down to very short wavelengths. It has al¬ 
ready been shown, in Eqn. (27), that these n—1 integrals 


are the n—1 power of a single integral J. The infinite ex¬ 
tent of the integration range, in turn, is a consequence of 
using continuum mechanics that has no intrinsic length 
scale, as opposed to real solids that have at least one mi¬ 
croscopic length scale given by the interatomic distance. 
In the problem we are considering in this paper, there is 
an additional length scale given by the length L of the 
dislocation segments, and the potential § has been ob¬ 
tained for wavelengths long compared to L. Accordingly, 
we now attempt replacing I-ji by the finite, but cut-off 
dependent, quantity /k(w,A) given by (43). How is A 


to be determined? We can not fix a precise value for the 
cutoff A, but we can provide and estimate by noticing 
that the divergence of I is controlled by the behaviour of 
G°{q) for |g| —>■ 00 . Now, G^{q) appears in the expression 
for the T-matrix because we use Eq. for the poten¬ 
tial. This is an approximate expression valid for wav es of 
wavelength longer than the dislocation segmenfpJiSI. Eor 
short wavelengths (large wave numbers) the dislocation 
dynamics can not be described by the exciting field in 
a reference position [for instance the middle of the dis¬ 
location segment as we take for long wavelength to get 
Eq.(|^]. Accordingly for large wavenumbers the poten¬ 
tial is different and that will change the behaviour of 
the integrand in I for |g| —> 00 . Therefore we chose an 
upper cut-off A = dwi that is proportional to the first 
resonance frequency wi = jV)^ with constant of 

proportionality d. Eq. (§ is expected to be valid up to 
these frequencies for a real material with dislocations of 
length L. 


Consider then (t°, given by (34), with I evaluated us¬ 
ing the cutoff A = dull, {0 < d < 1). Using (II) and 


(46) Arctanh^ « ^ we get 




aL 






with 


_ 8 ^352 

TT^ m 

oi = 1 -I- aLIoduJi 
d^ 

02 = 1 — oL/q—W i 


(48) 


(49) 

(50) 

(51) 


It is a straightforward exercise to verify that (48) leads 
to (471 in the limit ui ^ duJi ^ uji. 
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V. COMPARISON WITH THE 
RENORMALIZATION OF A DIRAC DELTA 
POTENTIAL IN A SCHRODINGER EQUATION 

The algebra of the problem we are dealing with has 
some similarities with the scattering of quantum (de 
Broglie) waves by an isotropic point obstacle as described 
by a Schrodinger equation in D dimensions: 


(v^ + fc^) V’(s’) = 


(52) 


This problem, because of its simplicity, has received some 
attention in the literature and it is instructive to compare 
it with our own. The interaction term in Eqn. (521 as 


our own Eqn. (§, is localized at a point. The main dif¬ 
ferences between this Schrodinger problem and our dis¬ 
location problem are: a) the former is scalar and the lat¬ 
ter vectorial and b) the potential in the dislocation case 
Eq.(|^ has two gradient operators. In the Schrodinger 
problem {k\V\k') = gj{2TT)^ i.e. a constant. In the 
dislocation problem {k\V\k') = —A{Mk- Mfc') as shown 
in Eq.(21). Note also that the two dimensional prob¬ 


lem of transverse (anti-plane) waves scattered by an in¬ 
finite straight screw dislocation is a scalar Schrodinger 
probleirP^ with a potential {k\V\k') proportional to k^. 
The analysis we review below also holds when the point 
obstacle is modelled by even derivatives of the delta 
functiorP^. 

When D = 1 the scattering problem associated to 
Eq. (521 poses no difficulty. When D = 2 the scat¬ 
tered cylindrical wave (and/or G°+) diverges logarith¬ 
mically at the origin thus creating a divergence in the 
T—matrix. Similarly, when D > 3 the scattered spheri¬ 
cal wave (and/or G®’*') diverges as l/r^“^ creating also 
a divergence in T. The idea that has been put forward to 
deal with this divergence is that the potential strength g 
in V(x) = gS^^^ {x) is not a measurable quantity but the 
scattering cross section (fc|T|fc) and possible bound state 
energies ar^^ In terms of the potential and the Green 
function one has 

{k\T\k') = {k\V\k') + {k\VG°V\k') + {k\VG°VG°V\k')-... 

a quantity we already computed for the dislocation prob¬ 
lem. Analogously, for the Schrodinger problem we find 




(53) 


with 


J'(fc) = I 


d^q 


(27r)^ — k'^ — ig 


Since J'(fc) diverges for I? > 2 one introduces a regulator 
A that cuts off the integral for jgl > A and the previous 
equation is interpreted as 


f 1 ^ 



I lim ( 

V27r 7 

/ A—>^00 \ 


/ 1 


J'{A,k) 


-1 


in which the unmeasurable parameter in the strength of 
the potential is turned into a function of A: g = g{A). 
Its form is then determined through the requirement that 
the scattering matrix be finite for some wave vector k = 
ko, because T{kQ) is an observable while g is not. Thus 

g{A) = g{T{ko),J'{A,ko))- 

To be explicit let us consider D = 5 

pA 


J'{A,k) = 


Q 5 


q^dq 


Iq q^ — k"^ — ig 

£5 A J 9 ^ 7 S A 7 ^ , n. 71 

+--1: Arc.anh-+.— 


(27r)5 
O.r; 


(54) 
fc^Tr' 


where Qs is the solid angle in D = 5 dimension. Note the 
similarity between the integral in J' with I from Eq.(38) 
in our own dislocation problem. 


Since it is not possible to cancel in Eq.(53) the differ¬ 


ent diverging functions of A that appear in Eq. (54) with 


1/g{A), because the former contains k dependent diverg¬ 
ing terms and the latter not, one finds from Eq.(53l that 
fimA_>.oo(fc|T|fc') = 0. In other words, there is no scatter¬ 
ing by a delta function potential when the interaction is 
governed by the Schrodinger equation. This is the case 
for all H > 4. The cases D = 2 and D = 3 are different. 
For instance with D = 3 


J'{A,k) = 


Q.q 


q^dq 


(27r)3 

1 

2'k‘^ 


) q^ — k"^ — ig 

A — fcArctanh^ -\- ik — 
A 2 


and if we measure the T-matrix at some fcp we get 
“At = [(27’')^T(fco)]"^ - lim J'(A,fco) 

K—^oo A—>-oo 

thus it is posible to define a renormalised coupling pa¬ 
rameter such that \/g{A) = I/^rn ~ 
divergence is absorbed in the unmeasurable coupling pa¬ 
rameter, which gives in Eq.( |53[ ) as A —)■ 00 a finite value 
for the T-matrix 


{k\T\k') = 


(^-^k^ 

V5RN 47r 


with (;rn a parameter to be determined by a measure¬ 
ment at some given k^. Similarly in T = 2 the diverging 
term of J' can be absorbed in the coupling parameter 
g{A). Moreover in the case of an attractive delta poten¬ 
tial 5 rn can be expressed in terms of the energy of the 
bound state (which is measured). In the case of repulsive 
delta potential the T-matrix (scattering cross section) 
must be measured at some value of k and the computed 
expression gives the value on the rest of the k axis. 

Let us now go back to the dislocation scattering prob¬ 
lem. Given the similarity between (431 and (54), one is 


tempted to perform the renormalisation procedure that 
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will lead to the conclusion that there is no scattering, as 
in the Schrodinger case for = 5. But the key difference 
is that the potential for the dislocation involves measur¬ 
able quantities such as elastic constants, frequency and 
Burgers vector. Thus we conclude that the finiteness of 
the cutoff value is unavoidable in this context. 

Let us finally remark that in the scattering problem of 
an anti-plane elastic wave by a screw dislocation in two 
dimensions the above analogy with the scattering of a 
Dirac delta holds, with D = A. 

VI. DISCUSSION 

The full Born series for the scattering amplitudes by a 
single edge dislocation segment as described by the po¬ 


tential Eq. is a geometric series that can be summed 
and expressed analytically. We have shown that apart 
form a multiplicative factor, the angular dependence of 
these scattering amplitudes are the same as in the lead¬ 
ing order Born approximation. The multiplicative factor 
is a frequency and cutoff dependent quantity, that di¬ 
verges if the cutoff goes to infinity. This divergence is 
an artifact of the perturbation series within a context of 
continuum mechanics, where there is no intrinsic, short- 
distance, length scale. At low frequencies, it has been 
shown that the seemingly divergent integral can be regu¬ 
larized to zero. As frequencies increase, the integral can 
be regularized through the introduction of a cut-off. fre¬ 
quency, proportional to the first resonant frequency of 
the dislocation segment. 


From an experimental point of view we can ask how the cutoff A can be determined. For this we note that from 
Eq.(43l in the large A limit (where ArcTanh(e) « e) we can get 


^) — In(^21 A) ^ I-ji{uj2,A) — I-ji , A) _ / uj^ — ~ <^2 \ 

w? - - wf A Vwl - wf wf - wl y ■ 

Thus, weighted differences of I at different frequencies can determine A. The consequence then appears to be that 
measuring the cross section at two different low frequencies will provide the appropriate value of A which is crucial 
to predict the cross section at any, third, frequency. And similarly for the mass operator. It appears then that in 
order to have a relation between frequency and wave vector for the coherent wave two parameters are needed, to be 
determined through measurement of uj{k) at two different points. 
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Appendix A: Second order Born approximation 

The second order Born approximation is given by 

Tj:^\k,k') = I dxd5?e-^^-^V,,{x)G%{x-x')Vik{x')e^^'-^' (Al) 

dQ ^in^mj^lp^kqkqknQmQpG 
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The second line is obtained using the expressions (j^ for the potential Vtj and (161 for the bare Green function 
The third line is obtained performing the integrations over x and x' with the aid of the delta functions. Finally, the 
last line uses the definition (23) of / and gives (221. 


Appendix B: Third order Born approximation 

The third order Born approximation is given by 
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Proceeding in the same way as for the second order gives Eqn. (251: 
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